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ABSTRACT 

(N ■ 

■ The Palomar Testbed Interferometer (PTI) is a long-baseline infrared interferometer 

located at Palomar Observatory, California. It was built as a testbed for interferometric 

{Sj ■ techniques applicable to the Keck Interferometer. First fringes were obtained in July 

1995. PTI implements a dual-star architecture, tracking two stars simultaneously for 

00 ! phase referencing and narrow-angle astrometry. The three fixed 40-cm apertures can 

be combined pair- wise to provide baselines to 110 m. The interferometer actively 
' tracks the white- light fringe using an array detector at 2.2 /im and active delay lines 

with a range of ±38 m. Laser metrology of the delay lines allows for servo control, 
and laser metrology of the complete optical path enables narrow-angle astrometric 
measurements. The instrument is highly automated, using a multiprocessing computer 
system for instrument control and sequencing. 

> ■ 

1. Introduction 

The Palomar Testbed Interferometer (PTI) is a long-baseline infrared interferometer installed 
at Palomar Observatory, California. It was developed by the Jet Propulsion Laboratory, California 
Institute of Technology, for NASA as a testbed for interferometric techniques applicable to 
the Keck Interferometer ( |Colavita et al. 1998 ). These include high-sensitivity direct-detection 



interferometry in the infrared with array detectors and dual-star interferometry with phase 
referencing for narrow-angle astrometry. However, the technology to implement PTI, including 
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active fringe tracking, optical delay lines, laser metrology, and real-time control, is applicable 
not only to the Keck Interferometer but also to other missions such as the Space Interferometry 
Mission (SIM) ( Unwin, Boden, Shao 1997| ). PTI also serves as a testbed for interferometric 
planning, operational, and data-processing and management tools applicable to both missions. 

Major development of PTI began in November 1992 with the commencement of funding from 
NASA under its TOPS program. The interferometer was installed at Palomar Observatory during 
the spring of 1995, and first fringes were obtained in July 1995. Subsequent work has focused on 
scientific investigations and engineering of the more advanced observational modes. Figure |l] is a 
photograph of the instrument. 

PTI builds on experience gained with the Mark III interferometer, which was a visible- 
wavelength interferometer installed at Mt. Wilson, California. The Mark III used active fringe 
tracking for wide-angle astrometry and parametric imaging (Shao et al. 1988; Shao & Colavita 
1992a). PTI differs from the Mark III in that it operates in the infrared, and in addition to 



parametric imaging supports narrow-angle interferometric astrometry implemented using phase 
referencing, described below. 



2. Narrow-angle interferometric astrometry 

The development of technology for detecting extrasolar planets is one of the motivations for 
building PTI. Most of the recent discoveries of extrasolar planets have been through measurement 
of radial velocities ( Mayor &: Qucloz 1995 ). This technique senses the doppler shift of the 



light from a star caused by the longitudinal component of an unseen companion's motion. 
A complementary technique to radial velocities is astrometry, which measures the wobble in 
the position of a star caused by the transverse component of a companion's motion. While 
high- accuracy global astrometry remains the purview of space missions, the detection of extrasolar 
planets through astrometry is primarily a narrow-angle problem, as the measurement can be made 
with respect to angularly-nearby reference stars. In this case the atmospheric limits to astrometric 
accuracy achievable with a long-baseline interferometer are several orders of magnitude below the 
wide-angle atmospheric limit ( [Shao fc Colavita 1992b ; |Colavita 1994 ). Figure || shows the increase 



in astrometric accuracy achievable in a differential measurement with long baselines and narrow 
fields. For interferometer baselines of ~100 m, and star separations of 15-20 arcsec, accuracies of 
tens of microarcseconds in an hour measurement are possible. 

To exploit this accuracy, both the parent star of the putative planet (henceforth primary 
star) as well as one or more astrometric reference stars (henceforth secondary stars) must be 
detectable by the interferometer. For the nearby sample, for which the astrometric technique is 
most sensitive, the primary star will be bright and provide an adequate signal-to-noise ratio within 
an TqTo coherence volume to track. However, angularly-nearby stars will generally be too faint to 
fringe track. But for secondary stars chosen within the isoplanatic patch of the primary star, the 
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primary star can be used as a phase reference (Shao & Colavita 1992a; |Quirrenbach et al. 1994| ) 



to synthesize a long coherent integration time, greatly increasing sensitivity. Phase referencing 
works best in the near-infrared, which provides a large isoplanatic patch and the ability to exploit 
large apertures. Combined with adaptive optics, which corrects wavefront aberrations across the 
individual pupils, the fringe stabilization provided by phase referencing allows for sensitivities 



similar to those of a space-based system ( Shao Colavita 1992a ). Sky coverage — the fraction of 



the sky over which a suitable phase reference is available — is a key metric. PTI, as a testbed with 
moderate apertures, has limited sky coverage. However, good sky coverage can be achieved with 
apertures in the 1.5-2.0-m range, as proposed for the Keck Interferometer outriggers ( |Colavita et 



al. 1998|) 



Combining phase referencing with simultaneous measurements of two stars is the basis of 
dual-star astrometry. Figure || shows the basic implementation. Light from separate stars in 
the field-of-view of each aperture is directed into separate interferometers sharing a common 
baseline. Internal path delays for each interferometer are measured with laser metrology to 
common fiducials at each aperture. These common fiducials "tie together" the two interferometers, 
allowing, for example, fringe-tracking errors on the primary star to be fed forward to the secondary 
star to implement phase referencing. The secondary interferometer, cophased by the primary 
interferometer, makes the actual astrometric measurement by switching between the secondary 
star and the primary star; the change in delay between the two stars is the astrometric observable. 

Because of the small field-of-view for the narrow-angle measurements, requirements on 
baseline knowledge are modest. The error 5 in a differential astrometric measurement over an 
angular separation 9 caused by a fractional baseline error of 5B /B is attenuated by the separation 
as 5 = 6(5B/B). For example, a 100 fim error in a 100 m baseline for astrometry over a 20 arcsec 
separation introduces errors of only ~20 fias. Thus, separate metrology of the baseline is generally 
unnecessary, and baseline solutions can be obtained from wide-angle astrometry using one of the 
interferometers QShao et al. 1990 ). 



3. Instrument overview 

PTI is located northeast of the 5-m Hale Telescope at Palomar Observatory, Palomar 
Mountain, California; Figure [l] is a photograph of the interferometer taken from the catwalk of the 
5-m telescope. The central building houses the optical delay lines, beam-combining laboratory, 
and the computer /control room. The three small shelters located north, south, and west of the 
central building contain the siderostats, telescopes, and dual-star feeds. Light is directed from the 
siderostat shelters through evacuated pipes to the central building. 

PTI is operated as a two-element interferometer, combining light from two of the three 
siderostats at one time. Initial observations combined light from the north and south piers for a 
110-m baseline which is oriented 20° east of north. The inclusion of the west pier allows additional 
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baselines. Table |T| summarizes some key instrument parameters. Figure |4| is a schematic of the 
beam train of PTI. Light at each aperture is collected by a 50-cm siderostat mirror, which directs 
the light to a fixed 40-cm telescope; a wide-angle camera adjacent to the telescope provides coarse 
star acquisition to get the stars within the field-of-view of the angle tracker in the beam-combining 
lab. The 40-cm telescope forms an image of the field, which is split into primary and secondary 
beams by the dual-star feed; fast-steering mirrors in the dual-star feed correct for wavefront tilt. 
The two beams from each aperture are directed through beam pipes to the central building where 
they are directed into a common long delay line, providing identical delays for the two beams. 
Local metrology of the delay line is used for servo control. Beam compressors reduce the beam 
diameters to a more convenient size before sending the primary star through a short delay line. 
This delay line introduces a differential delay between primary and secondary to enable phase 
referencing within the isoplanatic patch. PTI has separate beam combiners for the primary and 
secondary stars. In each the visible light (0.7-1.0 /im) is split off with a dichroic beamsplitter and 
directed to an angle-tracker which controls the fast-steering mirrors in the dual-star feeds. The 
infrared light (1.6-2.4 ^m) is combined in a Michelson beam combiner feeding one 128 x 128-pixel 
quadrant of a NICMOS-3 array. In dual-star mode, "constant-term" metrology (Sec. 11.2 monitors 
the entire optical path from beam combiner to dual-star feed. A distributed real-time control 
system orchestrates operation of the system. 



Design considerations Operation at 2.2 /xm is optimal for dual-star astrometry (Shao & 



Colavita 1992b| ), exploiting both the large coherence volume and increased isoplanatic patch 



available in the infrared. While the choice of a 40-cm aperture for PTI was established primarily 
on the basis of cost, it is also well matched to average atmospheric conditions. For nominal 
1 arcsec seeing, the coherence diameter ro at 2.2 /im is 53 cm; thus a d = 40-cm aperture is well 
corrected with only modest tip/tilt correction. 

The PTI design incorporates most of the usual system requirements. The layout is 
highly symmetric in order to minimize visibility losses due to differential image rotation and 
differential phase effects. To maintain high throughput, most optics use protected-silver coatings. 
Transmissive optics are generally antireflection-coated IR-transmissive fused silica, which work 
well out to 2.5 fim. Dust covers are used in most places to protect the optics, and a CC>2-snow 
gun is used to clean the siderostat mirrors. 

Atmospheric dispersion is relatively small over the K band (App. |C|), and thus air delay lines 
can be used, allowing for considerable simplicity of implementation. However, because of the 
long propagation paths in the delay lines and in beam transport, care must be taken to minimize 
turbulence, as discussed below. 

The system was designed to be highly modular. This allowed subsystem assembly and testing 
prior to installation at the site. For example, the pointing system was debugged on the roof of our 
lab at JPL, while closed-loop operation of the fringe tracker with the delay line was accomplished 
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in the basement of the same building. This modularity extended to the computer system, where 
subsystems were allocated their own real-time computer board, even if computational needs did 
not mandate it, in order to simplify the development process. 



4. Building and site 

The central building is a T-shaped modular insulated building with a length of 30.5 m. The 
top part of the building, 6.1 m x 18.3 m, houses the beam-combining lab and computer room; 
the latter is partitioned off from the rest of the building and located adjacent to the driveway in 
Fig. [l]. Figure || gives more detail. 

The beam-combining lab includes three optical tables. These optical tables are grouted onto 
0.45-m-diameter concrete posts on top of a central isolation pad. The center 1.2 m x 3.7 m 
table supports relay optics, beam compressors, metrology components, and the short delay lines. 
The other two 1.2 m x 2.4 m tables support 1.2 m x 2.4 m x 0.1 m breadboards on which 
are assembled the primary and secondary beam combiners; the use of breadboards simplified 
integration at Palomar, as most assembly could take place at JPL. 

The long part of the building, 4.9 m x 24.4 m, encloses the long delay lines. Two 22 m x 
0.6-m-deep I-beams provide the primary delay-line support. They attach via embedded J-bolts to 
the central isolation pad at the front and at 3 other isolated pads along their length. Box-beam 
crosspieces, 0.10 m x 0.15 m x 1.7 m, are located at 0.6 m intervals along the I-beam. Each delay 
line runs on a pair of 38.1-mm-diameter round steel rails with extruded aluminum supports which 
mount to the crosspieces; a center square rail is used for the friction drive. The support system 
allows for three delay lines, although only two are installed. 



A "building- in-a-building" , an approach first used on the SUSI interferometer ( Davis 1994| ), 



is implemented using 1.2 m x 2.4 m x 0.05 m rigid insulating panels surrounding the optical 
tables and long delay lines as shown in Fig. [|. The building-in-a-building allows the heating 
and air-conditioning units, which are mounted on isolated pads, to run continually, maintaining 
a constant interior temperature without affecting internal seeing. Additional thermal-control 
measures include the use of white paint on the building driveway and adjacent road; white-gravel 
is used in the vicinity of the siderostat shelters. 

The siderostat shelters house the siderostats, telescopes, and dual-star feeds. They are 
insulated buildings, 5.2 m x 2.7 m, elevated 1.5 m above the local grade. With the exception of 
one fixed wall, the roof and upper wall roll back on supports, as seen in the foreground shelter 
in Fig. |], leaving a 0.75-m-high lower wall. In operation, a windscreen made of rigid insulating 
panels surrounds the interior optics, and a separate cover encloses the dual-star feed. The shelter 
includes an electronics rack for the remote electronics, which is vented to the outside. 



The pipes that bring light from the siderostat shelters to the central building are evacuated 
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tubes capped with IR-transmissive fused silica windows. They terminate approximately 1-m 
from the siderostat shelters. During observations, when the shelter roofs are rolled back, short 
lengths of PVC pipe enclose the light paths from the end of the evacuated pipes to the edge of 
the dual-star feed cover to minimize local seeing. Similarly, at the termination of the evacuated 
pipes at the walls of the central building, short lengths of PVC pipe enclose the light path into the 
building-in-a-building. 

We note that the light pipes were originally simple non-evacuated insulated aluminum tubes. 
Sufficiently insulated, these tubes worked acceptably for propagating the 2-/im stellar beams. 



However, as part of a system upgrade in 1997 dWallace 1998j ), the tubes were replaced with 



evacuated pipes, which improved the propagation of the 633-nm constant-term metrology beams, 
which make a double-pass through the tubes. The driving disturbance in the air pipes was the 
temperature difference between the tops of the pipes, radiatively cooling to the sky, and the 
bottoms, which would be coupled to the warmer ground. While the insulation minimized the 
effect, some convection was seen at dusk when setting up the instrument; we suspect that the air 
pipes would be entirely acceptable if buried and away from diurnal disturbances. 



5. Siderostat system 

The siderostat mount is an alt-az gimbal moving a 50-cm zerodur flat. The mirror is 
supported on a 9-point wimetree suspension. Adjustable hard-points, opposed by spring plungers, 
locate the mirror surface with the altitude axis. The altitude bearing mounts can be shimmed 
to align the altitude and azimuth axes. The mount provides access to the back surface of the 
siderostat mirror through a hole in the mirror cell. A contacting depth gauge, encoded by a laser 
interferometer, was used to measure from a fixed point off the mount to the front or back of the 
siderostat mirror with the mount rotated to (az, el) coordinates of (0°, 0°), (180°, 0°), (0°, 180°), 
and (180°, 180°). Combined with the center thickness of the siderostat mirror, these measurements 
allow determination of the offset between the two axes and the mirror surface offset. (As the 
measurements are always made at near-normal incidence to the front or back surface of the mirror, 
tolerances on the precise mount angles at which measurements are made are relatively loose.) The 
mount is adjusted to keep deviations from an ideal siderostat to <100 /mi. Mirror edge location is 
via delrin tabs, and the unbalanced altitude load is compensated with a counterweight. 

The system is driven with a dc servo motor. It provides 0.027 arcsec motor steps and a 
maximum slew rate of ~l°/s. Direct axis encoding is provided by 6.48-arcsec optical encoders; the 
zero points of each axis are defined by magnetic proximity switches on the main wheel repeatable 
to less than one encoder step. Computer readable "soft" limit switches indicate extremes of travel 
while "hard" limits protect against software failures by disabling the siderostat drive. 
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Control aspects The real-time control system for the siderostats provides track and slew 
control based on sidereal targets from the instrument sequencer, as well as variety of other target 
types. Slew and track control use rate feedback to the axis indexers: large-angle slews use the 
encoders for closed-loop control, while tracking uses the indexer position; transition between 
modes is automatic based on the size of the tracking error. The controller accepts offsets from 
the acquisition cameras and angle trackers for star acquisition, and from the angle tracker for 
fast-steering- mirror desaturation. The translation from sidereal targets to mount position is via 
an 8-parameter mount model which incorporates the orientation of the axes, the zero points of 
the axes rotations, the orientation of the siderostat mirror in its cell, and the direction of the feed 
beam from the telescope. These parameters are determined via an off-line least-squares routine 
using pointing residuals to known targets; when properly calibrated, open-loop pointing is good 
to ~15 arcsec rms — well within the field-of-view of the acquisition system. The only nightly 
calibration required of the siderostat system is a homing process to establish the encoder zero 
points. Manual paddle controls are provided but are not required in normal operation. 

For observations in 1997 an accelerometer feedforward system was implemented to compensate 
for a small amount of piston attributable to the siderostats which was detectable in the fringe 
phase. The system uses an accelerometer^] mounted to the back of each siderostat mirror which 
senses acceleration along the mirror normal. Analog electronics filter and doubly integrate the 
accelerometer output to provide a measure of the mirror piston above 1.5 Hz. This signal is 
digitized at 2 kHz, scaled by 2 cos 6, where 9 is the angle between the siderostat normal and the 
telescope feedbeam, and is fed forward to the delay line for compensation. 



6. Acquisition system 

The siderostat mirror directs light to a 40-cm Cassegrain telescope. However, preceding the 
telescope, an elliptical mirror located in the shadow of the telescope secondary directs a fraction 
of the stellar beam to a wide-angle star acquisition system, as shown in Fig. ^. Given the long 
baselines and high magnification of the beam-transport system, accurate pointing is required in 
order for starlight to enter the field-of-view of the angle tracker inside the beam-combining lab. 

The acquisition system uses a commercial 12-bit thermoelectrically-cooled CCD camera^] 
which interfaces to the control system using GPIB. The camera is preceded by a 90-mm diameter, 
1300-mm focal-length telescope followed by a 2:1 teleconverter. These provide an angular scale 
of 1.8 arcsec/pixel over a field-of-view of 12x18 arcmin; smaller fields are selectable via subframe 
readout. 

The acquisition system incorporates a 70-mm corner cube for calibration. The corner cube 



5 Sunstrand QA3000 

6 Photometrics Star-1: format: 384x576; pixel size: 23 /jm; read noise: 20 e~; dark current: 15 e~/s at -45 C 
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is located at the edge of the elliptical mirror feeding the acquisition system so that a portion of 
a boresight laser injected in the beam-combining lab and fed backwards through the system is 
intercepted and directed into the acquisition camera. 

Control aspects The acquisition system implements a basic centroid operation to locate 
the star image in the field, using a 10x10 box about the brightest pixel. The exposure time 
is optimized if the intensity if out of defined bounds. Subframe readout is used after the first 
exposure to minimize image-transfer time. Feedback of centroids to the siderostats is used to 
automate the acquisition process. The acquisition system also provides information on apparent 
star magnitude to the angle tracker. 

The acquisition system provides several calibration modes. A primary boresight, done at the 
beginning of a night, records the image centroid of the primary boresight laser retroreflected off 
the acquisition corner cube; this location establishes the initial zero point for star centroiding. 
This zero point is updated between scans, once the angle tracker has acquired the star, using a 
star boresight, which records the current stellar position. The repeated star boresights remove 
drifts in the acquisition system as well as any systematic offsets due to aberrations in the primary 
laser image. For dual-star operations a secondary boresight measures the image centroid of the 
secondary boresight laser. This is used to servo the secondary star selector mirror in the dual-star 
feed to coalign the secondary and primary boresights. 



7. Collector telescope and dual-star feed 

The siderostat mirror feeds a 40-cm Cassegrain telescope. The telescope is mounted to a 
1.2 m x 2.4 m optical table, inclined 15° with respect to horizontal; the telescopes point 15° 
north of west at the siderostats ( |Colavita et al. 1994 ). The north and south telescopes employ 



a thermally-compensated truss ( |Krim 1975 ) with invar truss elements and aluminum end pieces, 



while the west telescope uses a rail-mounted secondary which is athermalized with an invar 
metering rod. The telescopes use an f/3 primary and provide an f/10 output beam at a back 
focal distance of 400 mm. The secondary diameter is 124 mm, with a secondary obscuration of 
140 mm. The wavefront quality of the telescopes is 30-40 nm rms. The 15-degree inclination of the 
telescopes improves the vignet-free field-of-regard of the siderostat. For sources at a declination of 
5 = 30°, the system is unvignetted for hour angles HA < 0; areal vignetting remains <10% for HA 
< 1.5 h and <20% for HA < 2.5 h. Figure || provides additional information on the field-of-regard 
of the system. 

The telescope feeds the dual-star feed (DSF), shown schematically in Fig. [?|. The DSF 
separates the light from the primary and secondary stars into separate beams and directs them 
to the beam-combining lab. The DSF also contains the fast-steering mirrors (FSMs) used to 
correct wavefront tilt. The DSF is constructed on a 0.9 m x 1.2 m optical breadboard mounted 
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horizontally behind the telescope. 

With reference to Fig. mirror Ml folds the inclined output beam of the telescope into the 
horizontal plane where it is focused onto the field separator M2 at a scale of 19.4 /xm/arcsec. 
The primary beam transmits through the field separator and is reflected by mirror M3 through 
a 25-mm hole in mirror M5 to parabola M4. The parabola has a focal length of 750 mm and 
recollimates the starlight into a 75-mm beam (compression of 5.33 with respect to the telescope 
primary). The collimated beam reflects off annular mirror M5 to fold mirror M6, which directs 
the primary beam toward the beam-combining lab. 

The large 75-mm collimated beam minimizes diffraction in the beam transport and also 
reduces delay-dependent phase and amplitude errors ([Tango fc Twiss 1974|) . The Rayleigh length 



d 2 /X for propagation at 2.2 fim is 2600 m; with ~100 m internal propagation distances, PTI 
operates at a small fraction of this value. 

While the primary beam transmits through the field separator M2, the secondary beam is 
reflected, and is recollimated and directed similarly to the primary beam with mirrors M4', M5', 
and M6'. Mirror M6' is referred to as the secondary star selector mirror (SSSM), and is located 
conjugate to the plane in the 40-cm beam between the siderostat and telescope which includes the 
telescope secondary. This plane defines the system input pupil, and also includes (physically or 
virtually) the constant-term fiducials, discussed in Sec. [LI. 2 . Tilt of the SSSM is used to select a 



secondary star from the telescope image on M2. An optional Lyot stop can be located adjacent to 
the SSSM. 

The SSSM is a flexure mount driven by two microstepped actuators; including the 5.33 
demagnification and a factor of 2 for reflection, its angular resolution on the sky is 0.015 arcsec. 
The SSSM is controlled by the siderostat control system, tracking the sidereal motion of the 
secondary star with respect to the primary star and accepting offsets from the acquisition and 
angle-tracker systems. 

Mirrors M4 and M4' are also flexure mounts and are driven with 15-/im PZT actuators for 
control of high-frequency tip/tilt sensed in the beam-combining lab. They have a range on the sky 
of ±10 arcsec. When tracking, only a small portion of this range is used. 

Several options are available for the field separator M2. For single-star observations, the field 
separator is just removed. For dual-star observations, and especially for engineering tests, an 
achromatic beamsplitter is convenient. For application in the converging beam we use a dielectric 
deposited on a 300-/im-thick substrate. For large magnitude differences between primary and 
secondary, the optimal field separator is a pinhole. However, to maintain metrology continuity, 
the pinhole must be deposited on a substrate incorporating a metrology beamsplitter coating. 
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8. Central optics and delay lines 

Figure |8| is a schematic of the central area of the beam-combining lab. The 75-mm collimated 
beams from the primary and secondary stars enter the central building as shown in Fig. ||. The 
light is directed by feed mirrors on the center table into the long delay lines, which occupy the 
long part of the central building. 



8.1. Long delay lines 

The long delay lines use a four-stage design nearly identical to one previously developed at 
JPL ( Colavita et al. 1991 ; Colavita, Hines, Sz Shao 1992; ). The optical assembly is an athermalized 



cat's eye implemented with a parabola and a flat. The parabola is 330x230 mm with a 1070-mm 
focal length. The flat at the system focus is mounted on a PZT for high-bandwidth servo control; 
a compensating PZT provides momentum balance. The optical assembly is suspended with flexure 
arms above an optics cart that rolls on steel wheels on the delay-line track. The optical assembly 
is driven with respect to the optics cart with a voice-coil motor for intermediate-level control. A 
separate motor cart contains a microstepper motor which provides coarse motion using a friction 
drive to the center rail; a second voice-coil motor drives the optics cart with respect to the motor 
cart. Control and drive signals are provided via a ribbon cable which runs in a cable tray from 
the front of the delay- line track, under the delay-line cart, around an idler pulley, and then to 
the motor cart. The idler pulley is on a third cart, which is tensioned via a cable from a torque 
motor at the back of the delay-line track. The delay line maintains delay jitter of 10-20 nm rms 
at speeds of tens of mm/sec. This is well in excess of PTI requirements of jitter < 50 nm for a 1% 
visibility reduction at 2.2 //m at a maximum sidereal rate of 5 mm/sec. 

PTI uses two long delay lines. The south ("active") delay line is of the type just described. 
The north ("passive") delay line includes just the motor stage and is moved only between 
observations. Thus quiet motion is not required and it can use a more rudimentary (and less 
expensive) optics cart that rides on ball bushings and has rigid connections to both the optical 
assembly and the motor cart. Residual motion of the passive delay line during an observation is 
sensed by laser metrology and servoed out by the active delay line. The delay-line tracks are 22-m 
long; allowing room for the various carts and some margin for deceleration from slew, the physical 
travel of each delay line is 19.15 m for a total optical delay using both delay lines of ±38.3 m. The 
delay-line travel was chosen subject to cost constraints, and sets the limits of sky coverage on the 
N-S baseline, shown graphically in Fig. ||. 

The delay lines incorporate a PZT for modulation which is mechanically in series with 
the servo PZT. This PZT is driven with a fringe-scanning waveform as required for fringe 
detection. The use of a separate PZT for modulation allows for simple distortion correction; a 
waveform-generator card provides a predistorted drive waveform to correct for PZT nonlinearities. 
If the modulation were implemented as an additional signal on the servo PZT, its nonlinearity 
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would be a function of the operating point of the servo. 

The long delay lines delay both the primary and secondary stars. With respect to the 
parabolic mirror, 75-mm primary and secondary beams enter at the upper left and upper right, 
and exit diagonally opposite the optical axis at the lower right and lower left. These delayed 
beams exit the delay lines and pass beneath the feed mirrors toward the beam compressors. 



8.2. Beam compressors 

The beam compressors reduce the starlight beam diameter by a factor of 3.75 from a transport 
diameter of 75 mm to 20 mm to reduce the size of the remaining back-end optics. They use an 
on-axis Gregorian design with confocal parabolas. The primary is 100 mm in diameter with a 
500-mm focal length. The compressor despace is maintained with invar rods for athermalization. 
The projected obscuration of the compressor is less than that of the main telescope secondary. 



8.3. Short delay lines 

After the beam compressors the primary beams pass through the short delay lines. These 
introduce a differential delay between the primary and secondary star. Thus, if L and I are the 
delays introduced by the long and short delay lines, the delays seen by the primary and secondary 
stars are L + I and L, respectively. With this delay-line configuration, pathlength modulation can 
be implemented on the long delay line, affecting both the primary and secondary stars, or on the 
short delay line, affecting only the primary star. 

The requirements on the small delay line are modest: for a ±30 arcsec range on the sky with a 
110-m baseline, the required delay is ±17 mm with a maximum sidereal rate of 1 /im/s. The south 
small delay line is active, with a PZT and motor, while the north one is completely passive and 
included for optical symmetry. The small active delay line implements a simple two-stage servo. 
It uses a PZT at the focus of an f/5 100- mm parabola and a microstepper-driven translation stage 
which moves the entire assembly. 



8.4. Beam combiner tables 

The 20-mm primary beams delayed by the short delay line are directed to the primary 
beam-combining table; the light from the secondary star is directed after its compressors directly 
to the secondary beam combiner. The beam-combiner tables include the optics and detectors for 



the angle tracker and fringe tracker. Figure |10| is a schematic of the primary beam combiner. The 
secondary table is similar, but uses a different fringe-detector configuration. The angle tracker 
uses light reflected from dichroic Bl into telescope LI, which feeds the angle-tracker detector. 
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The fringe tracker demodulates the interference fringe formed at beamsplitter B3. The angle 
tracker and fringe tracker are discussed in detail in Sees. ^ and [l0|, below. Additional components 
on the beam combiner table include a white-light source, which is collimated and injected into 
the combined beam via beamsplitter B5. The white-light source is used for boresighting and 
calibration of both the fringe tracker and angle tracker. A HeNe laser can be injected via 
beamsplitter B6 and is used for boresighting the acquisition system, as discussed in Sec. ||, as well 
as for manual boresighting and focusing. A microscope at the back of the system is useful for 
system alignment and test. Also indicated on the figure are the locations of the constant-term 
metrology injection and extraction optics (see Sec. 11.2). 



9. Angle tracker 

The angle tracker senses and controls the wavefront tilt of the interfering starlight beams. 
The system uses an angle sensor in the beam combiner that controls fast-steering mirrors in the 
DSF; identical systems are implemented on both the primary and secondary beam combiners. 
Angle sensing is over the bandpass 0.7-1.0 /im using a quad-cell detector implemented with 
photon-counting Si APD detectors; a chopping scheme is implemented which shares the angle 
sensor between the two input beams on each beam combiner. The bandpass is limited at long 
wavelengths by the bandgap of the Si detectors, and at short wavelengths by a long-pass filter 
which blocks scattered HeNe metrology light. 



Atmospheric issues Assuming a standard Kolmogorov turbulence model, the total wavefront 
variance across an aperture of diameter d is 1.03(d/Vo) 5//3 rads 2 ; with perfect tip/tilt correction, 
the variance is reduced to 0.134(d/r ) 5/3 ( |Noll 1976Q . If the angle tracker removes 90% of the 
tip/tilt variance, the residual wavefront error for one aperture is a = 0.37 rads rms. This results 
in a mean value for the squared fringe visibility of exp (-2a 2 ) = 0.76 (see [Tango fc Twiss 198CQ . 
To achieve this performance requires only modest angle-tracker bandwidths. For a wavefront tilt 
power spectrum H(f), and wind speed W, the product fH(f) peaks at 0. 075-0. 2bW/d depending 
on whether the wind is transverse or longitudinal ( Roddier et al. 1993| ); most of the tilt energy lies 
below these frequencies. Assuming W = 10 m/s and d = 40 cm, the peak is at 2-6 Hz. 



Description The angle tracker components are mounted on the beam combiner tables as shown 



in Fig. 10. The 20-mm collimated beams from each aperture enter from the left. Dichroic Bl 
reflects the bandpass from 0.7-1.0 /jm and transmits at 633 nm and 1.5-2.4 /im. The reflected 
light from Bl is directed to telescope LI (diameter = 90 mm; focal length (fl) = 1300 mm) which 
forms non-overlapping images of each beam on field lens L2 (fl = 38 mm). Mirrors Ml and Ml' 
are adjusted such that the beams cross at plane Rl prior to the field lens L2; this plane is imaged 
by L2 via fold mirror M2 onto mirror M3. Lens L3 (fl = 250 mm) relays the image from field lens 
L2 onto the lenslet array L4 which defines the 4 quadrants of the detector. L4 is portion of a 
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commercialQ replicated microlenslet array; the lenslet pitch is 400 fiioa with individual focal lengths 
of 3.2 mm. The lenslet is mounted to the front of a custom assembly containing 4 100-fim optical 
fibers at the focus of 4 of the lenslets; the fibers, with opaque jackets, feed the 4 detectors. Filter 
Fl is a long-pass filter which blocks scattered laser light at 633 nm. 

The detectors are commercial^] passive-quenched avalanche photodiode (APD) photon 
counters selected for low dark count (<5 counts/s). They provide a raw quantum efficiency 
of ~27% across 0.7-1.0 ^m and are usable to count rates of ~1 MHz. Both beams share the 
detectors. This is accomplished by mounting mirror M3 on a PZT tip/tilt actuator. M3 is driven 
by a 100-Hz square-wave signal, alternately placing images of each beam onto the lenslet array L4. 
The image scale at L4, referred to the sky, is 130 /xm/arcsec, and thus the lenslet array defines a 
quad-cell with a size of 6x6 arcsec. The separation of the chopped images is 8 arcsec. 

Control aspects The fundamental frame time of the system is 10 ms. Each frame includes 
a 4.3-ms integration on each beam; the dead time allows for settling of the chopper. Computer 
interface is through gated pulse counters. The fundamental track controller implements a 
conventional quad-cell algorithm and provides integral feedback to the FSMs in the DSF; control 
bandwidths to ~10 Hz can be implemented. Lower-bandwidth feedback is provided to the 
siderostat (or SSSM) for FSM desaturation. For faint sources, coadding of 10-ms frames can be 
selected. A spiral search algorithm, through direct handoff to the siderostat (or SSSM for the 
secondary angle tracker) is used for star acquisition. The search algorithm spirals out and then 
back to the origin using increasing search radii; the return to origin makes the algorithm robust 
against clouds and other beam blockages. The measured source intensity from the acquisition 
system is used to establish a detection threshold. 

The secondary angle tracker can also accept direct high-bandwidth tip/tilt feedforward from 
the primary angle tracker in order to implement an isoplanatic tilt correction from the primary 
star to the secondary one. This can be implemented in conjunction with low-bandwidth tracking 
on the secondary star itself to take out slow drifts. 

Calibration of the angle tracker uses two corner cubes to generate reference images. The 
white-light source injected at beamsplitter B5 propagates backwards through the system. It 
generates two beams which reflect off dichroic beamsplitter Bl, retroreflect from corner cubes C2 
and C2', and then transmit through Bl to produce images on the quad cells. The angle tracker 
adjusts the end-points of the chop waveform driving mirror M3 in order to center each image at 
the quad-cell center. This chop calibrate procedure is one of our initial nightly calibrations. 

The angle tracker also includes support for other less frequent calibrations, including 
measurement of detector quantum efficiency and calculation of the coordinate transformation 
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matrices between detector and FSM, detector and siderostat (or SSSM), and FSM and siderostat. 



10. Fringe tracker 

The fringe trackers on PTI measure the fringe parameters of the interfering light. The fringe 
phase is used to control the optical delay lines to track the atmospheric fringe motion. The fringe 
phase from the primary star can also be used to cophase the secondary star. More precisely, three 
modes are supported: 

single star In this case, the secondary channel is disabled. The primary fringe tracker works at 
high speed to follow the atmospheric motion. The short delay line is disabled, and control 
handoffs are given to the long delay line. While the nominal sample time is 10 ms, coadding 
adjacent samples can be used, under conditions of good seeing, to increase the signal-to-noise 
ratio. 

bright secondary This case adds fast tracking of the secondary fringe motion. Delay-line 

handoffs are trivially orthogonalized to allow independent tracking. As discussed in Sec. |8.3| , 
the long and short delay lines both contribute to the primary delay, while only the long delay 
line contributes to the secondary delay. Thus, if d p and d s are filtered delay updates for the 
primary and secondary, delay-line handoffs to the long and short delay lines are generated 
as 5L = d s and 51 = d p — d s . 

faint secondary In this phase-referenced mode, high-bandwidth control of the secondary fringe 
position comes from the primary fringe tracker, while low bandwidth control comes from the 
secondary star itself. Delay-line handoffs are orthogonalized as 5L = d p + d s and 51 = —d s . 



10.1. Atmospheric issues 

PTI uses coherent tracking of the white-light interference fringe, with lower bandwidth fringe 
centering obtained from the group delay. Continuity requirements applied to the modulo-27r 
fringe-phase measurements set the sample spacing. For 1-arcsec seeing, corresponding to a 
coherence diameter ro = 53 cm at 2.2 /im, and a constant wind speed W = 10 m/s, the fringe 
structure function is equal to one radian rms at coherence time ro,2 = 11 ms (see App. [A]). For 
a simple phase unwrapper, the sample spacing must be less than this value to minimize the 
probability of unwrapping errors, and we adopt a value of t =10 ms. With a more sophisticated 
phase unwrapper, as implemented on PTI, the effective coherence time increases to Tq 2 = 14 ms 
(see App. g). 

Fringe motion during the coherent integration time reduces fringe visibility squared as 
(Tango & Twiss 198C): exp (— T/Tq^) 5 / 3 - The quantity Tq^ is the coherence time associated with 
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fluctuations about the interval mean (see App. |A|), given by Tq^ = 0-815ro/W = 43 ms. For T 
= 6.75 ms (discussed below), the reduction in V 2 due to temporal effects is 5%. Finally, the 
two-aperture Greenwood frequency can be computed as fop. = (l-546r /Ty) _1 = 12 Hz. For a 
servo bandwidth f c = 10 Hz (1/10 of the 100-Hz sample rate), the residual fringe-tracker error is 
(/G, 2 // c ) 5/6 = 1.2 rads rms. 



10.2. Optical configuration 



The optical configuration of the primary beam combiner is shown in Fig. 1C. Dichroic 
beamsplitter Bl passes the infrared light from the two 20-mm input beams, and a flat fringe 
is formed at beamsplitter B3; compensator B2 accounts for the thickness of the beamsplitter 
substrate. The combined light from one side of the beamsplitter ( "white-light output" ) is directed 
via lens L5 (200 mm fl) and mirrors M5 and M6 to an external focus adjacent to mirror M10. 
The light from the other side of the beamsplitter ( "spectrometer output" ) is directed via dichroic 
beamsplitter B4 through a spatial filter assembly and external spectrometer to form a line focus 
on mirror M10 adjacent to the white-light spot. 

Dichroic beamsplitter B4 is 90% reflective from 1.5-2.4 //m; the residual transmission 
allows injection of the white-light and laser sources in the combined beam for calibration. The 
spatial-filter assembly is implemented using off-axis portions of parabolas M7 and M8 which focus 
the light into a single-mode fluoride-glass fibeij^] and recollimate the exiting light.0 As the spatial 
filter is implemented in collimated space with mirrors M7 and M8 on kinematic bases, spatial 
filtering can be removed if desired with only minor realignment of the system. To maximize 
tracking sensitivity, no explicit spatial filtering is used on the white-light channel, although some 
filtering results from the finite pixel size. 

Prism PI and lens L6 (200 mm fl) implement a low-resolution spectrometer. Prism PI can 
be chosen as desired, but is typically a simple 8-12° (physical) fused silica wedgeQ Lens L6 (like 
L5, L7, and L8) is a custom air-spaced achromat (fused silica + calcium fluoride), afocal from 
1.5-2.4 fj,m. 

The external focus at M10 is recollimated with lens L7 (100 mm fl) and relayed into a 
nitrogen dewar.p| The dewar contains a filter wheel in collimated space with K (2.00-2.40 ^m), 
K' (2.00-2.33 fj,m), and HK (1.5-2.40 /im) filters, plus a blank-off position. Pupil stops, nominally 



9 Le Verre Fluore: A c = 1.45 fim, NA = 0.176. 

10 The spatial filter was used for all observations in 1997, and several observations in 1996. 

11 The original spectrometer configuration used a single 16-deg fused silica prism; the configuration for 1997 reduced 
the dispersion with a -7.68-deg BK7 prism between M9 and L6; the text describes the current configuration. 

12 IR Labs ND8 
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10-mm diameter, are attached to the filters. A second 100-mm lens after the filter wheel relays 
light onto the infrared detector. 

The primary beam combiner, discussed so far, works at atmospheric rates. As it will typically 
be read-noise limited, the warm, slit-less spectrometer design described above is adequate. For 
the secondary beam combiner, phase referencing allows longer integration times and the infrared 
background becomes more important. Thus a cold spectrometer is implemented, which moves 



some of the functions shown in Fig. |10| into the secondary dewar. Figure |11| gives a schematic of 
the secondary dewar. Light from the white-light and spectrometer outputs of the beamsplitter are 
separately focused onto a slit behind the dewar window by 200-mm lenses LI and Ll', identical to 
those used on the primary beam combiner. An external in-line spatial filter assembly can be used 
on the spectrometer side as on the primary beam combiner. 

The entrance slit, nominally 40 fim wide, widens to 400 ^m for part of its length. The 
white-light beam enters through the wide part of the slit, is collimated by lens L3, passes through 
a filter wheel which includes the K, K', and HK filters and pupil stops, and is then relayed to 
the detector by lens L4. The spectrometer beam can use either the narrow or wide part of the 
slit. Because it enters off-axis, the spectrometer beam is intercepted by mirror Ml, is recollimated 
by lens L2, and dispersed by prism PI. Mirror M2 introduces a double pass, and is oriented to 
produce a line focus on the back of the slit adjacent to the white-light spot. This spectrum is then 
relayed through the filter wheel to the detector via lenses L3 and L4. For use on bright stars, the 
spectrometer prism PI can be removed, M2 realigned, and an external prism inserted prior to Ll'. 



10.3. Detector and electronics 

The detector for both beam combiners is an engineering-grade NICMOS-3 infrared array. 
One 128x128 quadrant is used. Pixels are 40x40 //m, and the external optics are configured so 
that the white-light spot and the spectrum are located adjacent to each other on a single line of 
the detector; the spectrometer resolution typically provides 5-8 pixels across the K band. 

Video amplification and clock buffering are provided with commercial electronics^ on pod 
boxes attached to the dewar. The video output drives a 16-bit, 500-kHz A/D converter which 
interfaces to the real-time system through a FIFO on a custom VME board. This board also 
includes a large downloadable memory that contains the clock patterns for the array. The clock 
patterns for real-time use are aperiodic, and are precomputed and downloaded during nightly 
system configuration. 

The NICMOS-3 detector provides a quantum efficiency of ~65% at 2.2 fj,m, and a read noise 
for a single 2-/xs correlated-double-sample read of ~32 e~. For use on PTI, multiple consecutive 
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2-/us reads are digitally averaged to decrease the effective read noise, typically to 11-16 e , as 
discussed below. 



10.4. Fringe demodulation and array clocking 

To implement fast tracking, fringe demodulation on PTI uses a fringe-scanning algorithm, 
similar to that employed on the Mark III Interferometer ( |5hao et al. 1988 ), but adopted for use 



with array detectors. The modulation waveform is a 100-Hz sawtooth pattern implemented on one 
of the delay lines. This unidirectional scan minimizes the number of reads required to measure the 
fringe parameters. The array timing is aperiodic, varying the integration time per pixel such that 
each pixel is scanned by exactly one wavelength. 

The timing for each 10-ms sample is identical. Each sample begins with a reset of the active 
and adjacent lines of the array followed by a short settling time. During this time the stroke 
waveform retraces from its value at the end of the previous sample. The first read begins typically 
1.5 ms into the sample and measures the reset pedestal on the white-light and spectrometer 
pixels. The stroke waveform begins its linear scan at this time, and successive reads occur after 
each quarter wave of modulation for each pixel. Each of these (nondestructive) reads is typically 
an average of 16-64 consecutive 2-fj,s subreads. With the dewar blanked off, the system at PTI 
achieves effective read noises for fringe demodulation of 24, 20, 16, 13, and 11 e~ at 2, 4, 8, 16, 
and 32 subreads. The non-white behavior most likely corresponds to 1/f noise in the detector. 

The array clocking involves some overhead associated with the array reset and settling time, 
although this does take place during the stroke retrace. Additional overhead comes from the 
shift register implementation on the NICMOS-3 chip, which does not have a clear, so that access 
to a previous pixel requires clocking out, reloading the shift register, and clocking back in.PI 
Fundamentally, matching the stroke to the wavelength means that the effective integration time at 
midband will be less than the longest integration time allowed. The present configuration achieves 
an integration time of 6.75 ms at 2.2 fim and 7.4 ms at 2.4 /im. 

With wider wavelength coverage to accommodate both H and K bands simultaneously, 
varying the bin widths with pixel becomes less practical. Appendix |D| discusses the corrections to 
the fringe quadratures necessary for use with fixed-width time bins. 



The PICNIC chip, which replaces NICMOS, has clearable shift registers and is resetable by row rather than by 
individual pixel. 
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10.5. Fringe acquisition and tracking 

The fast fringe-tracker algorithm is a variant of that used on the Mark I— III interferometers 
(Shao et al. 1988; Shao & Staelin 1980). Denote the five reads during each 10-ms sample 



time Cj , and di , where i = denotes the white- light pixel and i = 1 . . . R denote 

the R spectrometer pixels. The integrated flux in each quarter-wave time bin is calculated as 
Ai = cii — Zi, Bi = — a>i,Ci = Ci — bi, and Di = di — q . The raw fringe quadratures and total flux 
(in units of dn) are calculated from these values as 

Xi = Ai-d (1) 
Yi = Bi-Di (2) 
Ni = Ai + Bi + d + Di. (3) 

From these quantities we can estimate the fringe phase, visibility, and signal-to-noise ratio, but 
first it is necessary to correct for biases associated with the detection and readout process. 



A detailed discussion of the bias correction of these quantities is given elsewhere ( polavita 



1998). While explicit calibration measurements are interspersed with science observations, two 
initial calibrations are used to estimate the biases for use by the real-time system. The mean 
values B x , B , and B N on X, Y, and N are estimated from a low- level calibration with the 
system observing dark sky. Applying these bias corrections yields X = X — B x , Y = Y — B Y , 
and N = N — B N , from which an energy measure NUM is computed as 

NUM = A 2 + y 2 . (4) 

The squared quantity NUM contains biases attributable to read-noise and photon noise. The 
detector read-noise bias B rn is estimated from the value of NUM measured from a low-level 
calibration. The photon noise bias B pn takes the form B pn = kN, where k is the detector scale 
factor in dn/e~. The factor k is estimated from the increase in the bias on NUM with increasing 



light level. This high-level calibration uses the white-light source in Fig. fLOj, which is retroreflected 
back into the dewar with a corner cube placed at the exit of the primary table. 15 Thus NUM is 
corrected for these ac biases to yield 

NUM = NUM - B rn - B pn . (5) 

Using these bias-corrected quantities, instantaneous estimates of the fringe phase, fringe 
visibility, and fringe signal-to-noise ratio S (photon-noise limit) are computed as 

= tan- 1 ^ (6) 



15 The corner cubes CI, CI', while optically in the proper location, do not see much infrared light, as dichroic 
beamsplitter BI is nominally transmissive at that wavelength. 



= 2>™. (8) 
N 



With read noise, the signal-to-noise ratio (SNR) of the phase estimate is just 



4 N 2 V 2 

where a 2 ds is the effective double-correlated-sample read-noise variance. 

The fringe tracker supports three fundamental states: search, semilock, and lock. Transitions 
between states depend on the instantaneous value of the fringe signal-to-noise ratio S 2 , as 
computed above, or its average value S 2 , which is computed as the output of a sliding boxcar 
filter. A boxcar filter is selected for its finite impulse response (rather than a decaying-memory 
filter). Currently, a 15-point memory is used. 

In the search state the fringe tracker commands the delay line to implement a geometric 
spiral about the predicted fringe position. Typically, the search step is 4.4 fim every other sample, 
with an initial search out to +50 ^m, back through zero to -100 /im, etc. Transition from search 
to semilock occurs for S 2 > if. 

In the semilock state the fringe tracker attempts to track the white-light fringe. If S 2 > T 2 
after a time-out period, typically 100 ms, the fringe tracker enters the lock state; otherwise, it 
reverts back to the search state and continues the spiral pattern. 

In the lock state the fringe tracker tracks the white-light fringe. Loss of lock, and reversion to 
the search state and a restart of the search process about the last lock position, occurs if S 2 < T 2 . 
The transition thresholds T%, T2, and T3 are typically 6.0, 4.0, and 3.3; as implemented, they 
incorporate a weak dependence on photon flux to accommodate uncompensated biases at high 
light levels. 

The fringe-tracker servo is straightforward. The white-light phase at sample n, (f>(n), is 
unwrapped to produce <3?(n), which is filtered and used to update the delay-line position. While 
a fixed-bandwidth tracker is employed, a Kalman-filter phase unwrapper is used to improve 
performance. Let x{ n ) be the atmospheric phase process; a state-space model and estimator for 
X is presented in App. ^. The atmospheric phase is the sum of the unwrapped white-light phase 
plus the offset from the predicted sidereal position of the delay line, L, viz. xi n ) = &(n) + kL(n). 
The white-light phase 4>(n) is unwrapped about a prediction for the unwrapped phase at time n as 

$(n) = mod_ 7I - i7r [0(n) - <S(n|n - 1)] + %{n\n - 1), (10) 

where $(rt|n — 1) is computed with a Kalman-filter predictor of x( n ) as 

<fr(n|n — 1) = x(n|n — 1) — kL(n). (11) 
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The update for the Kalman filter uses the observation y{n) = <£(n) + kL(n). 

A simple integral fringe-tracker controller is implemented using the error signal $(n). 
Essentially, the delay-line track handoff is updated with a new rate at sample boundaries. The 
fringe-tracker gain is set to achieve a closed-loop bandwidth of 10 Hz. 



10.6. Fringe centering 

The K-band fringe envelope contains several fringes that have adequate SNR for tracking. 
For accurate phase referencing and astrometry, as well as for the best visibility measurements, it is 
necessary to choose the fringe closest to the center of the white- light fringe envelope, i.e., to work 
at zero group delay. While the group delay is readily estimated using the spectrometer channels, 
because of their lower SNR the group delay is not available at the full sample rate. Rather the 
group delay is estimated at a lower rate and used to correct for integral tracking errors. 

To improve their SNR, the spectrometer phasors Xi + jYi, i = 1 . . . R, are phase referenced 



to the white-light phase and coadded as described elsewhere ( Colavita 1998|) . The group delay 
is estimated from the peak bin of the complex Fourier transform of the phasors. Assume the 
spectrometer is linear in wavenumber from v\ = to vr = 1/Xr, for a synthetic wavelength 
A = (y\ — vr) -1 . With Ai = 2.4 /xm and A2 = 2.0 /im, A = 12 ^m. Let f p be the frequency of the 
peak in the power spectrum 



R 2 



j2nfi 



(12) 



The group delay g is then given as g = f p A(R — 1)/R. 



In addition to the group delay itself, a quality measure based on fluctuations in the power 
spectrum at other frequencies is also estimated. If the group delay is accepted, a track offset, 
typically limited to ±1 fringe, is applied to the delay-line target. Typical coherent integration 
times for estimating the group delay are 0.4-0.9 s. To simplify the data processing, track offsets 
are p ended to 0.5 s boundaries. 



10.7. Other fringe-tracker functions 

The previous discussion addressed the fast-tracking mode of the fringe tracker. The slow 
tracking mode will not be discussed in detail here, but a few comments are in order. With phase 
referencing and the cold secondary spectrometer, it should be possible to obtain background-limited 
detection. This requires a slower readout than described above to minimize the effect of read 
noise. At these low rates fringe scanning can be implemented using quarter-wave delay-line offsets 
rather than with continuous modulation. Another point is that cophasing quality is not limited to 
the tracker error, as set by the tracker bandwidth and Greenwood frequency. Essentially, as the 
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fringe tracker measures its tracking error accurately, that information can be used in an open-loop 
feedforward correction. The result is that cophasing quality is limited by time delays from the 
mean epoch of the fringe measurement to its application to the secondary channel. 

Besides tracking and acquisition, the fringe tracker implements initial and on-going 
calibrations to correct for biases in the fringe parameters as discussed above. Other functions 
include alignment modes for positioning the white-light and spectrometer images, and full-frame 
readout for testing and alignment. The fringe tracker also implements an FTS mode, which is 
useful for testing and for wavelength calibration of the spectrometer channels. 

10.8. Performance 

The photodetection probability for the fringe tracker on PTI is approximately 4%. This 
incorporates all effects, including the factor of 2 attributable to the use of one beamsplitter output 
for fringe tracking, the 6.75 ms / 10 ms duty cycle of the fringe detector; the spatial filtering 
attributable to the finite-size white-light pixel; detector quantum efficiency; and warm and cold 
transmissivity. For V 2 = 0.4 and an effective read noise of 12 e~ , an SNR for measurement of the 
fringe phase of 10 and 5 requires 1000 and 400 photons, respectively. The current fringe-tracking 
limiting magnitude is 4.5-5.0 K. 

11. Laser metrology 

PTI uses two types of laser metrology systems: local metrology of the long delay lines for 
closed-loop servo control, and end-to-end (constant term) metrology for astrometry and control of 
the short delay line. 

11.1. Long delay line metrology 

The long delay line metrology system uses a fairly conventional heterodyne system, similar 
to that used on the Mark III interferometer ( |Shao et al. 1988| ). The laser source plate uses a 
stabilized single-frequency HeNe laser and two acousto-optic modulators to generate a 2-MHz 
frequency offset between polarizations. The source plate feeds two metrology heads, one each for 
the active and passive long delay lines. For each metrology head, one polarization serves as a 
phase reference while the other polarization measures the delay-line position in double pass. The 
metrology light enters and exits through two additional apertures on the delay line front plate. 
Flat mirrors, rather than corner reflectors, are used for the metrology endpoints because of their 
tolerance to beam shear. A matched pair of optical wedges is used over the metrology apertures 
on the delay-line front plate so that the metrology light focuses to a slightly different location on 
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the delay-line secondary mirror; thus any laser scatter from secondary-mirror imperfections falls 
outside of the starlight field-of-view. 

The 2-MHz reference signal from the laser source plate and 2-MHz unknown signals from 
the active and passive long delay line metrology heads are digitized and fed to a custom laser 
fringe-counter card. The card counts integral and fractional fringes, providing 32 bits of integer 
and 8 bits of fraction, for a delay resolution with the double-pass implementation of 1.2 nm over a 
range of 1300 m. Limiting doppler excursions to ±1 MHz allows for delay slew rates of ±0.3 m/s. 



11.2. Constant-term metrology 

Assuming a reasonably stable optical train, the long delay line metrology described above is 
all that is needed for single-star amplitude observations. Drifts in the unmonitored optical paths, 
as well as drifts between the (slightly different) metrology and starlight paths through the delay 
line, are tracked out by the fringe servo and at worst affect only the acquisition time of the fringe 
tracker. However, for cophasing and astrometry, metrology of the entire starlight path, common 
mode with the starlight, is required. For wide-angle astrometry, both from the ground, and also 
from space flUnwin, Boden, &: Shao 1997 ), the interferometer baseline is defined by the location of 



a corner reflector located on the surface of the siderostat mirror. External metrology of the corner 
reflector monitors baseline changes, while internal (constant-term) metrology monitors internal 
optical paths from the beam combiners to the same point. However, as discussed in Sec. the 
requirements on baseline stability for narrow-angle astrometry are much less than for wide-angle 
astrometry, and it is possible to move the corner cube off the siderostat to a point in the collimated 
space between the siderostat and DSF. 

In the DSF design of Fig. 0, the SSSM is conjugate to the 40 cm beam from the siderostat 
at the plane containing the telescope secondary mirror. This pupil is thus a good location for a 
retrorefiector, as the metrology beam does not shear as the SSSM is adjusted. Alternately, this 
plane can be reimaged as shown in the figure. Lens LI in Fig. [7| forms an image of this plane (and 
the SSSM) in collimated space at the location of corner cube CC. A portion of the metrology light 
injected from the beam combiners along the primary and secondary paths thus retroreflects off of 
this common fiducial. 

The constant-term metrology system on PTI is a separate heterodyne system from the long 
delay line metrology system. The constant-term source plate uses a stabilized laser to generate 
two dual-polarization beams: one at a frequency difference of 110 kHz at a frequency offset of 
27 MHz from the laser frequency for the primary beam combiner, and a second with a frequency 
difference of 440 kHz at a frequency offset of —27 MHz from the laser frequency for the secondary 
combiner. The source plate samples each beam to generate 110 and 440 kHz reference signals. 



Initially, consider only the primary constant-term metrology. The 110-kHz dual- frequency 
beam is injected with a 5-mm mirror into the combined primary starlight beam in the space 
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between beamsplitters B4 and B5, as shown in Fig. 10 and retraces the starlight path of each 
arm to the DSF. Part of the light reflects off the star-separator beamsplitter to corner cube 
CC, retroreflects, and returns to the primary beam combiner where it is extracted by a dichroic 
beamsplitter located between beamsplitter B3 and lens L5. As the starlight beamsplitter is 
non-polarizing, it is necessary to add polarizers such that the s polarization traverses one arm and 
the p polarization the other. To achieve this, small polarizers are placed in each arm over the 
fraction of the pupil used by the metrology beam. Alternatively, polarizers could be used directly 
in front of the corner cubes in the DSF (although the polarization state has degraded somewhat 
by that point). 

The extracted metrology beam detected through a cross polarizer is a measure of the difference 
in the optical pathlengths between the two arms of the interferometer. With the exception of 
the thin polarizers, this implementation is fully common mode with the starlight, but it is not 
especially efficient with respect to laser light. The relatively low polarization difference frequency 
allows for more sensitive metrology electronics (than those used for the long delay lines) to ensure 
a good SNR. The constant-term metrology signal is digitized and measured with a digital fringe 
counter, similar to the one for the long delay line system. 

The secondary constant-term metrology system works identically to the primary one, using 
the same fiducials. However, examination of Fig. shows that when the SSSM is adjusted so 
that the primary and secondary beam combiners are boresighted (a calibration mode for dual-star 
astrometry), the primary metrology light leaks into the secondary path, and vice versa. The 
frequency arrangement above, with different difference and offset frequencies for primary and 
secondary, ensures that the only in-band leakage terms are at the frequency difference of the other 
metrology system. These leakage terms are filtered out of the analog signal with bandpass filters. 

For differential astrometry the primary and secondary constant-term measurements, CT p and 
CT S , are combined with the phase offsets from the primary and secondary fringe trackers to yield 
the astrometric delays 

d p = CT p + k~ l <f> p (13) 
d s = CTs + k- 1 ®^ (14) 

The difference of the two delays, A = d p — d s , is the narrow-angle astrometric observable. This 
difference is also used as the error signal for servo control of the small delay line. The difference A 
is calibrated by simultaneous measurements of the same (primary) star by both the primary and 
secondary fringe trackers. 



12. Real-time control system 



The real-time control system on PTI controls the hardware of the individual subsystems as well 
as orchestrating the interactions between subsystems and providing overall instrument sequencing. 
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A multiprocessing system is implemented to provide adequate computational resources. It also 
allowed a clean partitioning of functionality which enabled concurrent development. 

12.1. Hardware architecture 

Real-time control and sequencing on PTI is implemented with 7 VME-based 68040 single- 
board computers (SBCs) interfaced to a SUN workstation, as illustrated in Fig. |l2| ( pines 1994 ). 
All of the real-time computers use the Vx Works operating system and C/C++. Cardcage 1 
contains 4 SBCs, one each for control of the siderostat system, the angle-tracking system, the 
acquisition system, and overall instrument sequencing. While the angle-tracker detector interfaces 
directly to this cardcage using a custom photon-binning board, most of the controlled hardware 
is located at the siderostat piers. Small slave card cages at these locations are accessed via 
fiber-optic bus repeaters to allow low-latency hardware access. The remote cardcages contain the 
GPIB boards for interfacing to the acquisition system, indexers for controlling the siderostats and 
SSSMs, and D/A converters for driving the FSMs. 

The second cardcage contains a single SBC which controls the three servo-controlled delay 
lines: the long active, the long passive, and the short active (the short passive delay line has no 
actuators). This system implements a delay-line controller, described elsewhere ( |Colavita et al. 
[19911 |Colavita, Hines, k Shao 19921) . 

The third cardcage contains two SBCs, one each for the primary and secondary fringe 
trackers. The cardcage contains custom boards for clocking and interface to the IR arrays, as well 
as video cards for array display. 

Inter-cardcage communication is via a reflective memory (essentially an asynchronous shared 
memory). The software driver uses the reflective memory to implement FIFOs, double-buffers, 
and other structures for interprocess communication. Absolute time is available in each cardcage 
via an IRIG-B time-code reader; IRIG-B time comes from a GPS system. While each SBC uses 
ethernet for booting, only the sequencer implements a real-time ethernet connection to the SUN; 
all other subsystems interface to the sequencer (and each other) via the reflective memory. 

12.2. Software architecture 

The software architecture for each subsystem is similar, consisting of a number of periodically- 
scheduled tasks (1-2000 Hz) tied to GPS time, as well as event-driven tasks. Each subsystem 
implements a state machine characterized by modes which can be commanded, e.g., for the 
siderostat subsystem: track, home, etc., and associated states sequenced by the individual 
subsystem, e.g. track:slew, track:track, etc. Each subsystem provides a library of public functions 
that can be called by other subsystems. For the siderostat example, public functions are available 
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to pass pointing handoffs, provide tracking corrections, etc. These public functions use the 
interprocessor-communications layer to call private functions that provide the actual functionality. 
Common public functions for each subsystem support command interpretation, data recording, 
and status and error reporting. 

The sequencer is the only subsystem without explicit hardware. It provides high-level 
instrument sequencing and apparent-place calculations, automating observations to provide high 
observational throughput. The sequencer communicates with the hardware subsystems through 
their public functions, and implements a command relayer and data, status, and error servers. 
The sequencer interfaces with the SUN host using the RPC protocol over ethernet. RPCs are used 
to implement a command stream from the SUN to the sequencer, and to return error, status, and 
data streams. The SUN hosts a graphical user interface for instrument control and monitoring, as 
well as the data-recording system. 



12.3. User interface and sequencing 

The user interface provides several views into the system. The highest-level view is the 
observing window, which provides a summary of the status of each subsystem and control of 
the overall instrument sequencing. The instrument sequencer, a large state machine running on 
the sequencer SBC, provides fully automated observing (after nightly setup) based on star lists 
describing the targets and integration times. User control at this level is of the overall sequencer 
and data recorder. The sequencer is multithreaded, allowing concurrency to speed operations. For 
example, for automating a single-star observation, concurrent threads implement delay-line slew, 
and star acquisition and tracking; when the threads rendezvous, fringe acquisition can begin. 

For engineering and debugging, each subsystem provides a detailed window which displays 
a variety of status information and provides direct control of the subsystem mode, subsystem 
data recording, and access to many public functions, for example, to allow the passing of manual 
handoffs. All systems have a parameter structure that contains hardware mappings, servo 
parameters, etc., which can be updated from the engineering window. The default values of the 
parameter structures are updated from a free-form initialization file. 

The data-recording subsystem separates the combined data-recording stream, tags the data 
blocks, and generates separate files for each subsystem or instance of a subsystem. These files are 
recorded to disk and later transferred to CD-ROM for archiving. 
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A. Coherence time 

Coherence time can be defined in various ways (Buscher 1994). Let To,?, denote the 



structure-function definition of coherence time, viz. that sample spacing for which the phase 
difference between samples is one radian rms. The structure function depends on time as 
Di(t) = (t/ro,?) 5 / 3 . For i = 1, representing contributions from a single point on the wavefront (the 
usual adaptive-optics definition), to,i = 0.314ro/IF for coherence diameter tq and constant wind 
speed W ( [Buscher 1994] ) . For i = 2, applicable to interferometry, there are contributions from 
two apertures, and To, 2 = 0.207ro/W. Both definitions assume an outer scale much larger than 
ToiW, generally a good approximation. The definition of ro i2 assumes an interferometer baseline 
B S> TopW, valid for the long baselines of PTI. 

Let To t i to denote the variance definition of coherence time, viz. that time interval for which 
the phase fluctuations about the interval mean are one radian rms. It is given by Tq i = 1.235ro/W 
and To,2 = 0.815ro/VF, with time evolution erf = (T/To^) 5 / 3 . One- and two-aperture Greenwood 
frequencies ( |Greenwood 1977 ) can be computed similarly as fc,i = (2.343ro/VF) _1 and 



/g,2 = (1.546ro/W) . For a single-pole servo with bandwidth f c , the servo error variance is given 
by e 2 = (.A,,. ./;r :1 rad 2 . 



B. Estimating the fringe phase 

We model the time evolution of the total fringe phase x(t) by its structure function 
Dzif) = (t/To )2 ) 5 / 3 , valid for t <C B/W, where B is the interferometer baseline. The corresponding 
correlation function C{t) is given by (see [Ishimaru 197^ ) C 2 {t) = C 2 (0) - 0.5D 2 (*)- A state -space 



representation of the fringe phase process can be formed using Levinson's recursion to generate the 
autoregressive coefficients (see Kay fc Marple 198 1| ). In the limit of long baselines, C 2 (0) — > oo, 



and the one-state model for a sampling spacing t is given by 

Xn+l = X n + W n (Bl) 
Xn = X n , (B2) 

with plant noise = -D 2 (i)- For noiseless observations, the best one-step predictor of the fringe 
phase is Xn+i\n = Xn, i-e., just the previous point, with a prediction error given by the plant noise 
D 2 (t). 
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A two-state model is given by 
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with plant noise a 2 ^ = 0.655Z?2(t)- For noiseless observations, the best one-step predictor is 
Xn+\\n = Xn + 0.587(x n — Xn-i)-> with a prediction error given by the plant noise 0.655D 2 (t). 
Thus including a velocity term substantially reduces the prediction error. Incorporating the 
factor of 0.655 into the structure function, we can write the plant noise (and prediction error) as 
cr^ = D' 2 (t), where D' 2 (t) = (t/r^) 5 ^ 3 with an effective coherence time Tq 2 = 1.289ro,2- 



These matrix forms allow for a straightforward Kalman-filter implementation (see Gelb 1974 ) 
which can accommodate observation noise in measuring the fringe. Higher-order models provide 
little additional benefit: a three-state model has a prediction variance only 2% better than the 
two-state model. 



C. Atmospheric dispersion 

Use of an air delay line can introduce visibility losses, as the (dispersionless) vacuum delay of 
an off-zenith source is compensated by a dispersive air path. Fortunately, atmospheric dispersion 
is relatively small in the infrared. The second column of Table || gives the residual group delay as 
a function of wavelength when the phase delays are matched at 2.2 /im. These calculations assume 
a maximum air delay of 40 m at 0.8 atm. At A = 2.2 /im, the fringe envelope is offset 16 /im 
from the point of equal phase delays between the air and vacuum paths, with a spread of 12 /im 
between 2.0 and 2.4 /tm. Matching group delays between the air and vacuum paths at 2.2 fim 
yields the third column of Table which shows the residual phase as a function of wavelength 
about this point. Over the K band (2.0-2.4 /im), the quadratic phase residual yields a AV = 1.5% 
reduction in broadband (white- light) fringe visibility. 

For the narrow spectrometer channels, the spread in group delay across the band introduces 
only small visibility losses. For a spectrometer resolution R, the visibility loss attributable to an 
offset x from the peak of the fringe envelope for a particular channel is AV ~ (irx/(XR)) 2 /6; with 
a 5 fim offset and R = 40, the visibility reduction is ~0.5%. 



D. Fringe demodulation with fixed-width time bins 



When the effective pathlength-modulation stroke is not exactly equal to the wavelength, the 
fringe phasor estimates are distorted. However, they can be corrected in a straightforward manner. 
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Write the fringe pattern as 

f (</>) = - [X* cos <p + Y* sin 4> + N*}, (Dl) 

where s is the length of the pathlength modulation stroke, X* 2 + Y* 2 = (N*V) 2 , and 

<j) = taxT 1 (Y*/X*) - vr/2 (sin fringe). Define four time bins A = [-a/2, -a/4], B = [s/4,0], 

C = [0, s/4], and D = [a/4, a/2], and let X = A - C, Y = B - D, N = A + B + C + D. Then 

X = {l/s)[-a{s)X*-(3{s)Y*\ (D2) 

Y = {l/s)[+a(s)X* -0{s)Y*] (D3) 

N = (l/a)[+ 7 (a)X* + sN*], (D4) 

where 

a(a) = 2sin(s/4) - sin(a/2) (D5) 

/3(a) = l-cos(a/2) (D6) 

7(a) = 2sin(a/2). (D7) 

For a matched stroke (s = 2ir), a = (3 = 2. 

The measured values X, Y, and iV, can be corrected for the mismatched stroke as 

x COII = l+! + ^JL (D8) 

p a 

^corr = 7, (^^) 

p a 

N corr = N + j'(X-Y), (D10) 

where 7' = r y/(2a). While X corr , Y corr , and N COTT are unbiased, their variances are more complicated 
than for the case of a matched stroke, and there is generally some performance penalty. 
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Fig. 1. — Photograph of the Palomar Testbed Interferometer (PTI) at Palomar Observatory taken 
from the catwalk of the 5-m Hale telescope looking northeast. The driveway and adjacent road are 
painted white to minimize daytime heating. 

Fig. 2. — Differential astrometric accuracy vs. star separation in a one-hour integration for different 
baseline lengths (from Shao fc Colavita 1992b| ); atmospheric models providing 1/2- and 1.0-arcsec 



seeing are shown. These results assume an infinite outer scale, and better results are achieved 
when the baseline exceeds the outer scale, as would be expected with a 100-m baseline at most 
sites. Measurements with the Mark III interferometer of a 3.3-arcsec binary star flColavita 1994 ) 
are consistent with the model. 

Fig. 3. — Dual-star concept. 

Fig. 4. — PTI beam train. 

Fig. 5. — PTI building and overall layout. 

Fig. 6. — Schematic of the acquisition system. 

Fig. 7. — Schematic of the dual-star feed. Mirror M2 is the field separator; mirrors M4 and M4' 
are also fast-steering mirrors; mirror M6' is the secondary star selector mirror (SSSM). 

Fig. 8. — Detail of the beam-combining lab. Light paths for primary (solid line) and secondary 
(dotted line) stars are shown. The south delay line feed mirrors can be removed to allow injection 
of light from the west siderostat. The beams exiting the long delay lines pass beneath the delay 
line feed mirrors into the beam compressors. 

Fig. 9. — Sky coverage for PTI. The dashed lines bound the sky coverage for the N-S baseline as 
limited by the delay line range; the solid lines bound the sky coverage for the N-W baseline. The 
shaded areas on the figure indicate where areal vignetting of the beam exceeds 20%. 

Fig. 10. — Schematic of the primary beam-combiner table. 

Fig. 11. — Layout of the secondary dewar. 



Fig. 12 — 



PTI computer hardware architecture. 
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Table 1: PTI vital statistics 



Location 



Baselines 



Aperture 
Delay range 
Architecture 

Fringe-tracking wavelength 
Spectrometer wavelength 
Fringe sensor 

Angle-tracking wavelength 
Angle sensor 



Palomar Observatory, Palomar Mountain, CA 

(long: -116.8633°; lat: 33.3567°; elevation: 1687 m) 

NS: 110 m: ( 37.1 E, 103.3 N, -3.3 Z) 

NW: 86 m: (-81.7 E, -28.2 N, 3.1 Z) 

SW: 87 m: (-44.6 E, 75.1 N, -0.2 Z) 

40 cm 

± 38.3 m 

dual-star, active fringe tracking 
2.0-2.4 iim 
1.5-2.4 fim 

NICMOS-3 infrared array with fringe-scanning 

modulation, 10 ms frame rate 

0.7-1.0 fjtm 

Si APD quad cell 
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Table 2: Residual group and phase delays attributable to atmospheric dispersion at maximum 
delay-line position (see App. |C]). 

A (/xm) residual group delay (/im) residual phase (rads) 



1.6 


38.04 


11.05 


1.7 


32.74 


6.63 


1.8 


28.31 


3.70 


1.9 


24.56 


1.83 


2.0 


21.37 


0.72 


2.1 


18.62 


0.16 


2.2 


16.24 


0.00 


2.3 


14.16 


0.13 


2.4 


12.34 


0.47 
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